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ABSTRACT 

The goal of the current study is to identify scenarios for which thermal and moisture effects become 
significant in the loading of a composite structure. In the current work, a simple configuration was 
defined, and material properties were selected. A Fortran routine was created to automate the mesh 
generation process. The routine was used to create the models for the initial mesh refinement study. 
A combination of element length and width suitable for further studies was identified. Also, the 
effect of the overlap length and gap length on computed shear and through-thickness stresses along 
the bondline of the joints was studied for the mechanical load case. Further, the influence of 
neighboring gaps and overlaps on these joint stresses was studied and was found to be negligible. 
The results suggest that for an initial study it is sufficient to focus on one configuration with fixed 
overlap and gap lengths to study the effects of mechanical, thermal and moisture loading and 
combinations thereof on computed joint stresses. 


NOMENCLATURE 


Eu 

G a 

Gi c 

Guc 

h 


t 

T, 

To 

To, dry 
To, wet 

a i 

Pi 

8 


Yl 


Length of resin gap in joint 
Young’s modulus 
Shear modulus 
Mode I fracture toughness 
Mode II fracture toughness 
Ply thickness 

Model/specimen thickness 
Glass transition temperature 
Stress free temperature 
Stress free temperature, dry 
Stress free temperature, wet 
Coefficient of thermal expansion 

Coefficient of moisture expansion (swelling coeff.) 
Length of overlap in joint 
Shear strains 


R. Krueger, National Institute of Aerospace, 100 Exploration Way, Hampton, VA, 23666. This work was performed at 
the Durability, Damage Tolerance and Reliability Branch, MS 188E, NASA Langley Research Center, Hampton, VA, 
23681, USA 


l 



e ii 

Normal strains 

*7 

Exponent of B-K failure criterion 

V H 

Poisson’s ratio 

P 

Density (specific weight) 

°u 

Normal stresses 

T.. 

Shear stresses 


BACKGROUND 

Environmental factors that affect the performance of composite laminates include temperature, 
moisture content and hygrothermal history. Early studies showed that residual stress relaxation due 
to hygrothermal history altered the static mechanical strain necessary to cause matrix cracking in 
experiments [1]. Depending on material type and hygrothermal history, this effect was shown to 
either increase or decrease the resistance to matrix cracking. Matrix crack growth in laminates 
subjected to creep load and severe environments may require the development of viscoelastic failure 
criteria [1]. The viscoelastic response of laminated composites under the influence of mechanical 
and hygrothermal loads was studied analytically [2], Application of the analysis included the effect 
of residual thermal stress relaxation upon creep in laminated composite plates under mechanically 
applied loads [3]. Generalized plane-strain finite elements were also formulated using linear 
viscoelastic theory. Analyses of changes in interlaminar stresses at the free edge of multidirectional 
laminates subjected to cyclic load conditions indicated that several variables can affect peak 
amplitude in interlaminar stresses and cycle-to-cycle variations in stress state. These influencing 
parameters were determined to be frequency, stress ratio, relative dwell time and environment [4], 


PROBLEM OVERVIEW 

An overview of a scenario where combined mechanical and thermal load effects become significant 
was discussed in detail in reference 1 and is presented in Figure 1. For the pure elastic case, the 
transverse residual stress in the ply is tensile and the stresses increase as the temperature drops from 
the stress free temperature To, dry to an ambient value (see dashed red line along path © ->• (D). 
During moisture absorption, the transverse residual stress drops and may become compression (see 
dashed red line along path (2)). An increase in temperature will increase the compressive stress (see 
dashed black line along path (D). The stress increases during moisture desorption (see dashed black 
line along path ®) and the cycle starts from the beginning (see dashed red line along path ©). A 
stress relaxation, however, may occur caused by accelerated viscoelastic effects due to an increase 
in temperature and moisture content (see solid red line along path (3)’). As the temperature drops, the 
tensile transverse residual stress continues to increase beyond the original residual tensile stress (see 
solid red line path along (E>) [1], For the viscoelastic case, therefore, the presence of the thermal load 
could make the stresses worse due to the effect of the previous creep. 

These combined thermal and moisture effects may become significant in the loading of a composite 
fuselage structure where several individual components are connected with longitudinal joints, as 
shown in Figure 2. Since airframe manufactures have recently become interested in joint 
configurations, which may be designed with multiple overlaps and resin gaps as shown in Figure 3, 
it has become important to understand the stresses at the bondline. Overlap length, 8, and resin gap 
size, c, shown in Figure 3a may also vary. As the overlap length, 8, approaches zero (8^ 0) the 
design resembles the traditional splice joint, as shown in Figure 3b. 
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The current study, therefore aims at identifying scenarios where thermal and moisture effects 
become significant in the loading of a composite structure. Linear and geometrically non-linear 
finite element analyses of simple and more complex examples are required first, to study the single 
and combined effects of mechanical, thermal and hygral loading on the computed strains and 
stresses. The results from this initial study of simple specimens and load scenarios will be used as 
guidance to define the next steps. These next steps may include additional analysis of more complex 
specimens under different load combinations and more advanced viscoelastic analysis, as well as 
testing. 

The focus of the present work was to study the effect of overlap length and gap length on shear and 
through-thickness stresses along the bondline of staggered overlap joints using linear and 
geometrically non-linear finite element analyses. First, the material properties obtained from the 
literature are discussed. Second, the selected specimen configuration is introduced. Third, the finite 
element analysis is described, including a Fortran routine that was created to quickly generate two- 
and three-dimensional models of different configurations of interest. Fourth, the results of an initial 
mesh refinement study are discussed for the mechanical load case. Fifth, the effects of the overlap 
length and gap length on computed shear and through-thickness stresses at the bondline are 
examined. A summary and outlook for future steps are provided. The description of the Fortran 
routine and the required user input, as well as an Abaqus/Standard® input deck, are given in the 
appendix. 


MATERIAL DATA 

For the current study, material data for carbon/epoxy IM7/8552 were used since it is a commonly 
used material in the aerospace industry. Elastic constants, cure temperature and hygrothermal 
expansion coefficients for this material were collected from data available in the open literature [5- 
20]. Details about the data gathered are discussed in appendix A. A recent discussion indicated that 
the stress free temperature may actually be higher than published T g values [21, 22], However, this 
may be process dependent and may require a separate investigation and separate material 
characterization. Additional material characterization may also become necessary to generate 
constants required as input data for viscoelastic finite element analysis. 


SELECTED CONFIGURATION 

A simple configuration to study the effects of overlap (S) and gap (c) lengths on overlap stresses 
is shown in Figure 4. This configuration represents the idealized local section of an overlap 
bounded by adjacent plies representative of the center section of a specimen. It is based on prior 
experience with similar analyses [21, 22], For this simple configuration, a layup of [45, 90, -45, 
0]s was considered, where the center 0° plies were overlapped [23]. An overlap length, <5, varying 
between 6.35 mm and 63.5 mm and a gap size, c, varying between 1.27 mm and 12.7 mm were 
proposed [23]. These dimensions are of interest to the aerospace industry based on manufacturing 
considerations for advanced tooling and hand-layup tolerances [22-23]. For modeling purposes, 
boundary conditions were considered such that the section is constrained at one end (w=v=0) and 
displacements (w=const.) are applied at the other, thereby approximating tensile loading commonly 
applied to such joints. This model was used first to perform a mesh refinement study to determine 
the degree of refinement adequate to capture the stress gradients in the overlap region as discussed 
later. Second, an investigation of the effects of changing the dimensions <5 and c in the range 
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mentioned above were investigated. Two-dimensional plane-strain models were used first for linear 
and geometrically non-linear finite element analysis. 

A more complex configuration with 16 plies [45, 90, -45, 0]2 S and several staggered overlaps and 
gaps is shown in Figure 5. This configuration may be investigated at a later point including three- 
dimensional finite element analysis as well as specimen manufacturing and testing. Details of the 
specimen configuration such as layup, length and width would be defined after a successful study of 
the simple configuration shown in Figure 4. 


ANALYSIS 

Finite element modeling 

The commercial finite element software Abaqus/Standard® [24] was used to analyze the mechanical 
load cases and is capable of handling the thermal, moisture and combined load cases as well. For the 
current study only two-dimensional analyses were performed using higher order plane-stress CPS8 
solid elements to model the specimen. 

For more advanced analysis, a quasi-static procedure in Abaqus/Standard® may be used to analyze 
viscoelastic material and two-layer viscoplastic material behavior. The quasi-static procedure can 
also be used to analyze quasi-static creep and swelling problems, which occur over fairly long time 
periods (rate-dependent plasticity: creep and swelling). 

Fortran routine for automatic mesh generation 

A FORTRAN routine was developed that will automatically create an Abaqus/Standard® [24] input 
deck (.inp file) and thus simplify the variation in dimensions and mesh size for specimen types of 
interest. Specimen dimensions and material data including ply orientation and stacking sequence are 
provided by the user in a simple data file that is read by the routine at the beginning of the 
execution. Gaps are modeled by using the resin material properties in the respective regions of the 
finite element mesh. Examples of two-dimensional models of specimen types that can be created are 
shown in Figure 6. The mesh density in the thickness direction and locally in the length direction 
can easily be changed, which allows for mesh refinement and convergence studies. Different 
geometries of interest can also be generated such as a single lap-shear specimen shown in Figure 6b 
- 6d or a double lap-shear specimen as shown in Figure 6e. By changing the Abaqus/Standard® 
element type and specifying additional input parameters for the third dimension, the meshes shown 
in Figures 7 can easily be created. An example mesh for the single lap-shear specimen is shown in 
Figure 7a. In addition to changing the mesh density in the thickness direction and locally in the 
length direction, the routine also allows mesh refinement across the width of the specimen as shown 
in Figure 7b for the example of a double lap-shear specimen. The desired input for mechanical loads 
and boundary conditions can also be generated automatically as shown in Figure 8. An example of 
the prescribed displacements and fixed boundary conditions is shown in Figure 8a. Appropriate 
Abaqus/Standard® input for the thermal and moisture distributions, which are represented as fields, 
can also be generated. Currently only constant fields over the entire domain, as shown in Figures 8b 
and c, can be automatically generated. The details of the routine and a brief guide on its use are 
provided in appendix B. An example for user input is also provided together with a 
Abaqus/Standard® input file that was created using the routine. 
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Generalized plane-strain approach 

Several assumptions are possible when using two-dimensional finite element models. A plane-stress 
model enforces the out-of-plane stresses to be zero (crzz=Txz= T yz =0) and allows the displacement to 
be the free parameter. The plane-strain model, on the other hand, enforces the out-of-plane strains to 
be zero (£ Z z=7xz= 7yz=0), which excessively constrains the plies. The effect is most marked for 45° 
plies because of their high in-plane Poisson’s ratio, while it is small for 0° and 90° plies. Hence, the 
two dimensional plane-stress and plane-strain conditions may serve as upper and lower limits 
compared to the three-dimensional solution [25]. 

A different approach to account for three-dimensional width effects was suggested in reference 26 
such that modified material properties are entered in a two-dimensional model to create a 
generalized plane-strain case. The generalized plane-strain case is regarded as an intermediate state 
between plane-strain and plane-stress, and assumes that the out-of-plane normal strain £ 77 =-v\_e xx , 
where vl is the laminate Poisson’s ratio, and the shear strain y xz = y yz =0. With these assumptions, ply 
stiffnesses were calculated for each ply angle in the specimen and entered as material properties for 
the two-dimensional plane-stress models. Calculated properties are listed in Table II. 

Application of load and boundary conditions 

The load and boundary conditions shown in Figures 4 and 8 were used to simulate a region that is 
part of a larger section but would also mimic the condition of a test specimen clamped in hydraulic 
grips with applied constant displacement. The displacement was applied to the model such that the 
specimen was subjected to an axial tensile strain £**=0.01. It was planned to apply a constant 
temperature over the entire specimen (as shown in Figure 8b) and use the difference between cure 
and room temperature (AT=157 K) as input. For the moisture loading, a constant moisture 
distribution over the specimen was assumed (as shown in Figure 8c) to avoid the calculation of 
moisture diffusion through the thickness at this point. A literature study suggests that 1% increase in 
mass or weight due to moisture absorption is a reasonable assumption to simulate a saturated 
condition in the specimen [27-28]. 

Current finite element model 

The configuration selected earlier (shown in Figure 4a) was used for mesh refinement studies. For 
this first set of analyses, an overlap length of <5= 25.4 mm and a gap length of c=2.54 mm were 
chosen. Abaqus/Standard® plane-stress CPS8 solid elements were used to model a section that was 
3d=76.2 mm long, as shown in Figure 9a. A local displacement «=0.762 mm (0.03”) was applied at 
the left edge of the model to simulate an axial strain £^=0.01 in the section, modeled as shown in 
Figure 9b. The right edges of the section were fixed as shown in Figure 9c. The detail of the area of 
interest is shown in Figure 10 where the gap is modeled as a region c=2.54 mm (0.1”) with resin 
properties. The element length, Ax, shown in Figure 10b was varied during the refinement study. 
The paths defined in Abaqus/Standard® along which the stresses were computed (results shown 
later) were also included in Figure 10b (light blue lines). Due to a change in material properties, the 
computed stresses differ from one side of the path to the other when crossing the interface between 
plies. Therefore, stresses had to be computed for both sides of an interface for the same path, as 
shown in Figures 10c and d. 

The mesh refinement study was performed to determine the mesh required for subsequent analysis. 
Meshes were refined locally in the length and thickness directions as shown in Figure 11. A 
selection of local meshes used in the area of interest are shown in Figure 11. The element length Ax 
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was varied between Ax=0. 847mm for the initial coarse mesh, shown in Figure lib, and 
Ac=0. 026mm for the most refined mesh, shown in Figure lie. Simultaneously, the element height 
Ay was varied between Ay=0. 1 5mm (one ply thickness) for the initial coarse mesh shown in Figure 
1 lb and Ay=0.025mm for the most refined mesh, shown in Figure lie. 

The element length and height combinations modeled were: 

• Ar=0.847mm, Ay=0. 15mm - Figure 1 lb 

• Ar=0.423mm, Ay=0.075mm - Figure 1 lc 

• Ar=0.212mm, 4y=0.05mm 

• Ac=0. 1 06mm, Ay=0.0375mm - Figure 1 Id 

• Ar=0.053mm, Ay=0.03mm 

• Ac=0.026mm, Ay=0.025mm - Figure lie 

Based on the results of the mesh refinement study (discussed in the next section) element 
dimensions Ar=0. 026mm and Ay=0.025mm where chosen for the remainder of the study. To 
determine the effects of overlap length, 3, and a gap length, c, three additional configurations were 
studied: 


• <5=25.4 mm, c=1.27 mm 

• (5=25.4 mm, c=12.7 mm 

• (5=6.35 mm, c=2.54 mm 

Another model was created to simulate the case of multiple overlaps and gaps. The overlap length 
of (5=25.4 mm and a gap length of c=2.54 were held fixed but two additional overlaps and gaps were 
added to the configuration as show in Figure 12 to study if the presence of neighboring overlaps and 
gaps had a major effect on the results calculated for the initial 0° overlap in the center. A minimal 
influence of neighboring overlaps and gaps would indicate that, for a preliminary study, it is 
sufficient to study the effects of mechanical, thermal and moisture load on a simple configuration 
with one overlap, as shown in Figure 4a. 


ANALYSIS RESULTS 
Mesh refinement study 

Computed deformation 

The deformed finite element model is shown in Figure 13 for an applied displacement i/=0.762 mm 
which was used to simulate an axial tensile strain £ xv =0.0 1 . The deformation of the entire section is 
shown in Figure 13a where the undeformed outline is depicted as a dashed line. The local 
deformations shown in Figures 13b to 13e correspond to the refined meshes shown previously in 
Figures 1 lb to 1 le. The results suggest that a refinement through the thickness may also be required 
to capture the local deformation in the -45° ply above the resin region. 

Computed stresses 

Stresses were plotted along the pre-defined paths shown in Figure 10. Due to a change in material 
properties, the computed stresses differ on the sides of the path when crossing the interface between 
plies. Therefore, stresses were always plotted for both sides of an interface for the same path. For 
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the current study, the focus is on the shear stresses x xv and the through-thickness stresses a yv plotted 
on both interfaces along path 2. Path 2 was selected because the axial load, which is carried in the 
individual 0° plies, is transferred via shear along the overlap and therefore high stresses are expected 
in this region. However, further investigation is necessary to determine whether the stress along 
paths 1 and 3 could become critical under certain load conditions [21, 22], 

The computed shear stresses r xy and through-thickness stresses o yy along path 2 are plotted in 
Figures 14 through 17. The results indicate that refined meshes are required to fully capture the 
steep stress gradient in the vicinity of the overlap. The zone of computed high stresses and gradients 
is small and extends to about 2.5 to 3 mm on either side of the juncture of the resin zone and the 
terminated 0° ply (x=2.54 mm). Farther away from this location, coarser meshes yield the same 
results as fine meshes. Shear stresses x xy quickly drop to zero as expected, as shown in Figures 14 
and 16. The through-thickness stresses o yy reach a constant value outside this region, as shown in 
Figures 15 and 17. Based on these results, element dimensions Ar=0. 026mm and 4F=0. 025mm 
were chosen for the remainder of the study for the resin gap and the zone of peak stresses. Larger 
elements were used in the area where constant stresses were observed. For the current study, only 
trends were monitored and computed stress values were not compared to any failure stress as 
suggested in references 21 and 22. 

Computed stresses for linear and geometric non-linear analysis 

A comparison between stresses computed from geometrically non-linear analysis and results 
obtained from linear analysis are shown in Figures 18 through 21. Except for the through-thickness 
stresses, o m along the upper interface of path 2 (shown in Figure 19), the results follow the same 
trend and are in excellent agreement. The through-thickness stresses <% in the top 0° ply computed 
from the geometrically non-linear analysis decrease close to the 0° ply/resin boundary (green 
crosses) while the results from linear analysis sharply increase (red Xs). Away from the 0° ply/resin 
boundary, both results are in excellent agreement with each other. This discrepancy requires further 
investigation. However, based on the excellent agreement at all other locations, it was decided to 
use linear analysis for the remainder of the study since it allows simple superposition of results. 

Effects of overlap and gap lengths on computed bondline stresses 

To determine the effects of overlap length, 6, and gap length, c, additional configurations were 
analyzed. The individual results for computed shear stresses r xy and through-thickness stresses o w 
along path 2 are plotted in Figures 22 through 37. The left edge of the resin pocket was chosen as 
the x=0 location for all plots as above. The results confirm that to the right of the 0° ply/resin 
boundary and to the left of the 45° ply/resin boundary (x<0) the computed shear stresses x xv decrease 
to zero and the through-thickness stresses <%, reach a constant value over a short distance, x. The 
study therefore focused on the results at the zone near the 0° ply/resin boundary where the 
computed stresses have the highest values, as suggested in references 21 and 22. 

A comparison between stresses from all configurations is shown in Figures 38 through 41. The 0° 
ply/resin boundary was chosen as the x=0 location for the plots so that all stress profiles could be 
combined in one graph for comparison. The computed shear stresses r xy shown in Figures 38 and 40 
are in excellent agreement for all configurations considered. Varying the overlap length, <5, and gap 
length, c, did not affect the results. The computed through-thickness stresses o yy plotted in 
Figures 39 and 41 exhibit some sensitivity, in particular for the small values of the overlap length, 
8. For the short overlap length, some transverse deflection of the specimen was observed which may 
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cause the difference in computed through-thickness stresses o yy . This difference requires further 
investigation. Nevertheless, the results suggest that for an initial study, it is sufficient to focus on 
one configuration with a fixed overlap length, 6, and a gap length, c, to study the effects of 
mechanical, thermal and moisture loading and combinations thereof which is in agreement with 
previous findings [21, 22], 


SUMMARY 

The objective of the current study was to identify scenarios in which thermal and moisture effects 
become significant in the loading of a composite structure. In the current work, a simple 
configuration was defined, and material properties were selected. A Fortran routine was created to 
automate the mesh generation process. The routine was used to create the models for the initial 
mesh refinement study. A combination of element length and width suitable for further studies was 
identified. Also the effect of the overlap length and gap length on computed stresses at the overlap 
interface and resin gap was studied for the mechanical load case. Further, the influence of 
neighboring gaps and overlaps on the results was studied and was found to be negligible. The results 
suggest that for an initial study, it is sufficient to focus on one configuration with a fixed overlap 
length and a gap length to study the effects of mechanical, thermal and moisture loading and 
combinations thereof on overlap stresses. A user guide and example input values, as well as a listing 
of an Abaqus/Standard® input file, are provided in the appendix. 


FUTURE STEPS 

With a configuration selected, the mesh refinement study complete, and a critical area identified, the 
study now needs to focus on identifying scenarios where thermal and moisture effects become 
significant in the loading of a composite structure. Linear and geometrically non-linear finite 
element analyses need to be carried out to study the effects of 

1. Mechanical load only 

2. Thermal load only 

3. Moisture only 

4. Combination of mechanical load + moisture 

5. Combination of mechanical load + thermal load 

6. Combination of thermal load + moisture 

7. Combination of all three load cases 

Depending on the results of the initial analyses, more advanced analyses may be required. These 
analyses may include 

• Viscoelastic analysis to determine where the stress can be additive 

• Time dependent material properties 

• More complex example geometries 

• Specific structural geometries and load histories that have been identified to be a potential 
problem. 

The specimen configurations to be manufactured and tested, as well as the test matrix, should be 
decided upon after the successful completion of the proposed analyses. 
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TABLE I. MATERIAL PROPERTIES USED FOR CURRENT STUDY 


IM7/8552 Unidirectional Graphite/Epoxy Prepreg 
fiber volume approximately 60%, ply thickness h= 0.15 mm 

En = 161 GPa 

E22 =11.38 GPa 

£33 = 11.38 GPa 

vi 2 = 0.32 

V13 = 0.32 

y23 = 0.45 

G\2 = 5.2 GPa 

G 13 = 5.2 GPa 

G 23 = 3.9 GPa 

Glass transition temperature 

Tz,dry= 177 °C 

Tg,wet =143 °C 


Coefficients of thermal expansion 

a =- 9 10 7 1/°C 

a 2 =2.88 10' 5 1/°C 

a, =2.88 10' 5 1/°C 

Coefficients of moisture expansion 

l\ = 0.0 1/%M 

/3 2 = 3310 10 _ 6 1/%M 



TABLE II. MATERIAL PROPERTIES USED FOR 2D FE ANALYSES 


IM7/8552 Unidirectional Graphite/Epoxy Prepreg, ply thickness h=0A5 mm 
Input properties calculated using Minguef s generalized plane-strain approach [27] 

0 ° - ply 

£ll = 161 GPa 
vi 2 = 0.467 

E22 = 14.39 GPa 

G 12 = 5.2 GPa 

90° - ply 

£ll =9.81 GPa 
vi 2 = 0.461 

E22 = 12.42 GPa 

G 12 = 3.9 GPa 

±45° - ply 

£ll =32.4 GPa 
vi 2 = 0.464 

£22 = 13.8 GPa 

Gi 2 = 4.55 GPa 
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transvers ply residual stress 



© Initial transverse residual stress is tensile 

© Transverse residual stress increases as temperature drops to an ambient value 
© Transverse residual stress drops due to moisture absorption and becomes compressive 
© Compressive stress increases as temperatures increases 

©' Stress relaxation caused by accelerated viscoselastic effects due to increase in temperature and 
moisture content 
@ Desorption of moisture 

© A drop in temperature results in an increase in original tensile residual stress 


Figure 1 .A schematic diagram of viscoelastic changes in the stress free temperature due to 

environmental history [1]. 
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Figure 2. Fuselage components with longitudial joints. 




VZZZZZ 2 ZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZ1 



5: overlap 
c: resin gap 


a. Staggering of overlaps . 


ZZZZZZZZZZZZZZZZZZZZZZZZZZZL 


1 


zzzzzzzzzzzzzzzzzzzzzz^. 




ezzzzzzzzzzzzzzzzzz^^^^^^^^ 

^ •2223ZZZZZZZZZZZZZZZZZZZZZZZZZZ2 . . 


as 5-*-0, splice joint 


b. Splice joint. 


Figure 3. Joint configurations with overlaps and gaps. 
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IM7/8552 
h=0.15 mm 


y-v 


x,u 


u=const. 




yyyyyyyyyyyyyyyyyyyyyyyyyyyy 


IZZZZZZZ2\L 


V77777777777777777^JJJ3Jyysi. 


dimensions not to scale 


6.35 mm < 5 < 63.5 mm 
1.27 mm < c < 12.7 mm 
t= 1.2 mm 


A ■07“ 

/ — e 7 8 

',-2\ 


19 mm < 38 < 190 mm 


v.~e 3 


09 


e 2 


i 

u=v=0 


6^ 45° 

0 2 = 90° 

0 3 =-45° 

0 4 = 0°/gap/-45 c 

0 5 =-45°/gap/O° 

0 6 =-45° 

0 7 = 90° 

0«= 45° 


Figure 4. Simple configuration to study the effects of overlap and gap lengths 

and peform mesh refinement. 
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Figure 5. More complex configuration of staggered overlaps and gaps. 
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x,u 

a. Mesh for initial study of configurations in Figures 4 and 5. 


y-v 



x,u 



i 


y.v 



x,u 


c. Single lap-shear specimen with increased overlap. 



x,u 

d. Single lap- shear specimen with increased arm thickness. 



Figure 6. Two-dimensional finite element models created by FORTRAN routine. 
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a. Boundary conditions and mechanical load. 


y.v 



TEMP 

(Avg: 75%) 

+4.500e+02 
+4.352e+02 
+4.205e+02 
+4.058e+02 
+3,910e+02 
+3. 762e+02 
+3.615e+02 
+3.468e+02 
+3.320e+02 
+3. 172e+02 
+3.025e+02 
+2,878e+02 
+2.730e+02 



y.v 



b. Thermal load case. 


FV1 

(Avg: 75%) 

B +1.000S-02 
+ 1.0000-02 
+9.167e-03 
+8.333e-03 
+7.500e-03 
+6.667e-03 
+5.833e-03 
+5.G00e-03 
+4,167e-03 
+3.333e- 03 
+2.500e-03 
+ 1.667e-03 
+8.333e-04 
+0.000e+00 



c. Moisture load case. 


Figure 8. Load, boundary conditions and field input created by FORTRAN routine. 
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x,u 

a. Section modeled. 



b. Local detail of left edge and boundary condition. 



u=v=0 


) 



















| 




y> v 

I 

i 

(i 

P 



) 


£ 




c. Local detail of right edge and boundary condition. 


Figure 9. Finite element model. 
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y.v 


38=76. 2mm (3.0") 


a. Section modeled. 



b. Local detail of resin region c and ply overlap 5. 



























resin 0° 



c. Path 2 on upper interface. 



Figure 10. Definition of paths for stress plots . 
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y.v 


y- v 


y.v 


y.v 


y.v 


x,u 


a. Section modeled. 


x,u 



b. Local detail of coarse mesh bx=0. 847mm, Ay =0.15mm. 


x,u 



c. Local detail of coarse mesh bx=0. 423mm, Ay =0. 075mm. 


x,u 



d. Local detail of coarse mesh Ax=0. 106mm, by =0.037 5mm. 


x,u 



e. Local detail of coarse mesh bx=0. 026mm, Ay =0. 025mm. 
Figure 1 1 . Finite element meshes used for mesh refinement study. 
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Figure 12. Specimen configuration with multiple overlaps. 



u=0. 762mm (0.03") undeformed outline 


a. Deformed section modeled. 


Y.v 


x,u 


b. Local detail of deformed model for coarse mesh Ax=0. 847mm. 



c. Local detail of deformed model for Ax =0.42 3 mm. 



d. Local detail of deformed model for Ax=0. 106mm. 


y.v 


x,u 


e. Local detail of deformed model for refined mesh Ax=0. 026mm. 
Figure 13. Deformed finite element model. 
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Figure 14. Stress x along path 2 in upper 0° ply interface for 6=25. 4mm, c=2.54mm. 
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Figure 15. Stress o along path 2 in upper 0° ply interface for 6=25. 4mm, c=2.54mm. 
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Figure 16. Stress x along path 2 in lower 0° ply interface for 8=25. 4mm, c=2.54mm. 
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Figure 17. Stress o along path 2 in lower 0° ply interface for 8=2 5. 4mm, c=2.54mm. 
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Figure 18. Stress x along path 2 in upper 0° ply interface for 8=25. 4mm, c=2.54mm. 
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Figure 19. Stress o along path 2 in upper 0° ply interface for 5=2 5. 4mm, c=2.54mm. 
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Figure 20. Stress x along path 2 in lower 0° ply interface for 8=25. 4mm, c=2.54mm. 
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Figure 21. Stress o along path 2 in lower 0° ply interface for 8=2 5. 4mm, c=2.54mm. 
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Figure 22. Stress x along path 2 in upper 0° ply interface for 6=25. 4mm, c=1.27mm. 
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Figure 23. Stress o along path 2 in upper 0° ply interface for 6=25. 4mm, c=1.27mm. 
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Figure 24. Stress x along path 2 in lower 0° ply interface for 8=25. 4mm, c= 1.27mm. 
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Figure 25. Stress o along path 2 in lower 0° ply interface for 8=25. 4mm, c= 1.27mm. 
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Figure 26. Stress x along path 2 in upper 0° ply interface for 8=25. 4mm, c=12. 7mm. 
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Figure 27. Stress o along path 2 in upper 0° ply interface for 8=2 5. 4mm, c=12. 7mm. 

yy 


29 


stress t , 

xy 

N/mm 2 



distance along path x, mm 

Figure 28. Stress x along path 2 in lower 0° ply interface for 8=25. 4mm, c=12. 7mm. 
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Figure 29. Stress o along path 2 in lower 0° ply interface for 8=2 5. 4mm, c=12. 7mm. 

yy 


30 


500 



-2.0 -1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 


distance along path x, mm 

Figure 30. Stress x along path 2 in upper 0° ply interface for 6=6. 35mm, c=2.54mm. 
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Figure 31. Stress o along path 2 in upper 0° ply interface for 6=6. 3 5mm, c=2.54mm. 
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Figure 32. Stress x along path 2 in lower 0° ply interface for 8=6. 35mm, c=2.54mm. 
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Figure 33. Stress o along path 2 in lower 0° ply interface for 8=6. 35mm, c=2.54mm. 
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Figure 34. Stress r along path 2 in upper 0° ply interface for mutiple overlaps. 
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Figure 35. Stress o along path 2 in upper 0° ply interface for multiple overlaps. 
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Figure 36. Stress x along path 2 in lower 0° ply interface for multiple overlaps. 
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Figure 37. Stress o along path 2 in lower 0° ply interface for multiple overlaps. 
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Figure 38. Stress r along path 2 in upper 0° ply interface for all configurations studied. 
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Figure 39. Stress o along path 2 in upper 0° ply interface for all configurations studied. 
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Figure 40. Stress r along path 2 in lower 0° ply interface for all configurations studied. 
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Figure 4 1 . Stress o along path 2 in lower 0° ply interface for all configurations studied. 
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APPENDIX A: SELECTION OF MATERIAL DATA 

A study was performed to gather material data published in the open literature. The goal was to 
obtain a complete set of input data for one material which could be used in a finite element analysis. 
The material set should include the ply thickness, h, the elastic moduli, E. p Gk, and the glass 
transition temperature or stress free temperature, T g , as well as the hygrothermal coefficients of 
expansion cx r and fi t . None of the sources queried [5-22] provided a full set of data for one material. 
The most complete data sets were found for IM7/8552 listed in Table A-I. The data in the last 
column (boldface blue font) were selected for the current study. The pooled data are listed in the 
right column of Table A-I (boldface blue) and are identical to the values shown in Table I. Data for 
other materials are listed in Tables A-II through A-IV for completeness. The coefficients for 
moisture expansion, />(, for IM7/8552 could not be found in the open literature. Therefore, the data 
for IM7/8551-7 listed in Table A-IV were used together with the T g for that material (highlighted in 
boldface blue) and included in the right column of Table A-I to make up a complete set of input 
data. 


TABLE A-I. MATERIAL PROPERTIES FOR IM7/8552 FROM DIFFERENT SOURCES 


Property 

[5] 

[6] 

[7] 

[8] 

[9,11] 

[9,11] 

[10] 

[12, 13] 

Selected 

data 

Fiber Vol. % 



59.1 

60 

57.4 

60.3 

57.7 

57.3 


h [mm] 

0.15 


0.125 

0.125 

0.185 

0.152 

0.131 

0.183 

0.15 

Eli [GPa] 

150 

161 

171.4 

165 


162 

164 

162 

161 

Ell [GPa] 

10.7 

11.38 

9.08 

9.00 



12 

8.97 

11.38 

E33 [GPa] 

10.7 

11.38 


9.00 

9.3 




11.38 

V12 

0.33 

0.32 

0.32 

0.34 


0.311 


0.316 

0.32 

V13 

0.33 

0.32 


0.34 





0.32 

V23 

0.45 

0.45 


0.5 

0.475 




0.45 

Gu [GPa] 

4.8 

5.2 

5.3 

5.6 

4.4 

5 


4.7 

5.2 

Gu [GPa] 

4.8 

5.2 


5.6 





5.2 

G 13 [GPa] 

3.4 

3.9 


2.8 

2.3 




3.9 

Gi c [kJ/m 2 ] 


0.212 

0.277 

0.277 





0.212 

G// c [kJ/m 2 ] 


0.774 

0.788 

0.788 





0.774 

n 


2.1 

1.634 






2.1 

Eg (dry) [°C] 



AT=155 

AT=154 

177 

200 

200 

208 

177 ++ 

Tg (wet) [°C] 





143 

154 

154 

161 

143 

a, [1/°C] 

-9 10 7 


-5.5 lO -6 

-1 KT 6 





-9 10 7 

a 2 [l/°C] 

2.88 10 5 


2.5 10' 5 

1.8 1(T 5 





2.88 10 5 

a 3 [1/°C] 

2.88 10 5 



1.8 1(T 5 





2.88 10 5 

[1/%M] 









0 

P 2 [1/%M] 









3310 10 -6 

p [g/cm 3 ] 

1.55 




1.58 

1.58 

1.57 

1.58 

1.55 

[ 2 /m 2 ] 





190 

160 


190 













*The material properties are given with reference to the ply coordinate axes where index 1 1 denotes the ply principal 
axis that coincides with the direction of maximum in-plane Young’s modulus (fiber direction). Index 22 denotes the 
direction transverse to the fiber in the plane of the lamina and index 33 the direction perpendicular to the plane of the 
lamina. 
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TABLE A-II. MATERIAL PROPERTIES FROM DIFFERENT SOURCES 


Property 

[14] 

[15] 

AS4/8552 

[16, 17] 

grph/epxy 

[1] 

IM6/3501-6 

[3] T300/F- 
155 

Fiber Vol. % 




57 

57 

h [mm] 



0.127 



Eli [GPa] 



135 

149.2 

127.3 

Ell [GPa] 



11 

8.8 

7.9 

E33 [GPa] 






V12 



0.301 

0.34 

0.29 

V13 






v 23 




0.42 

0.42 

Gn [GPa] 



5.8 

4.5 

3.5 

Gi3 [GPa] 






d3 [GPa] 




3.1 

2.8 

Gj c [kJ/m 2 ] 






G// c [kJ/m 2 ] 






» 1 






Eg (dry) [°C] 


215 

177^ 

177^ 
171. Hf 

121^ 

119.3H 

Eg (wet) [°C] 






a, [1/°C] 

0 

6 10' 7 

-0.41 10 -6 

-0.54 10 -6 

0.36 10 -6 

a 2 [l/°C] 

3 lO' 5 

2.86 10 -5 

26.8 10 -6 

30.6 10 -6 

27.9 10 -6 

«,[1/°C] 

3 10' 5 





/?, [1/%M] 



0 

0 

0 

P 2 [1/%M] 



5560 10 -6 

2860 10' 6 

2530 10' 6 

P [g/cm 3 ] 






[g/m 2 ] 













*The material properties are given with reference to the ply coordinate axes where index 1 1 denotes the ply principal 
axis that coincides with the direction of maximum in-plane Young’s modulus (fiber direction). Index 22 denotes the 
direction transverse to the fiber in the plane of the lamina and index 33 the direction perpendicular to the plane of the 
lamina. 

+ To conform to experimental observations, it is assumed there is no expansion below 0.4% absorbed moisture [15] 

++ cure temperature 
^stress-free temperature 
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TABLE A-III. MATERIAL PROPERTIES FROM DIFFERENT SOURCES 


Property 

[1] 

IM7/8551-7 

[3] 

GY70/339 

[3] 

T300/934 

Fiber Vol. % 

57 



h [mm] 




En [GPa] 

151.6 

290 

145 

Ell [GPa] 

7.9 

6.1 

10.3 

£33 [GPa] 




V12 

0.34 

0.31 

0.29 

V13 




V23 

0.42 



G\1 [GPa] 

3.9 

4.1 

4.8 

G 13 [GPa] 




G23 [GPa] 

2.6 



Gic [kJ/m 2 ] 




G// c [kJ/m 2 ] 




n 




Tg (dry) [°C] 

177^ 

151.71 



Tg (wet) [°C] 




a, [1/°C] 

0.36 10 -6 

-0.72 1 (T 6 

-0.09 10 -6 

a 2 [l/°C] 

28.8 10 -6 

25.2 10 -6 

29.7 10 -6 

«,[1/°C] 




P t [1/%M] 

0 



P 2 [1/%M] 

3310 10 -6 

+ 5400 1 0 6 

+ 5400 1 0 6 

P [g/cm 3 ] 




[g/m 2 ] 









*The material properties are given with reference to the ply coordinate axes where index 1 1 denotes the ply principal 
axis that coincides with the direction of maximum in-plane Young’s modulus (fiber direction). Index 22 denotes the 
direction transverse to the fiber in the plane of the lamina and index 33 the direction perpendicular to the plane of the 
lamina. 

+ To conform to experimental observations, it is assumed there is no expansion below 0.4% absorbed moisture [15] 

++ cure temperature 
^stress-free temperature 
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TABLE A-IV. MATERIAL PROPERTIES FROM DIFFERENT SOURCES 


Property 

[18, 19] 

T300/5208 

[18, 19] 

AS4/3501 

[20] 

IM7/8551-7 

Fiber Vol. % 




h [mm] 



0.127 

Eli [GPa] 



165.5 

Ell [GPa] 



8.3 

E33 [GPa] 




V12 



0.18 

V13 



9.3 10' 3 

v 23 




Gu [GPa] 



6.2 

Gi3 [GPa] 




G 13 [GPa] 




Gj c [kJ/m 2 ] 




G// c [kJ/m 2 ] 




r\ 




Eg (dry) [°C] 


111 


Eg (wet) [°C] 




a, [1/K] 

0.02 10 -6 

-0.3 10 -6 

-0.04 10 -6 

a 2 [ 1/K] 

22.5 10 -6 

28.1 10 -6 

18 10 -6 

a 3 [l/K] 




0, [1/M] 

0.0 

0.0 

0.01 

[1/M] 

0.6 

0.44 

0.35 

P [g/cm 3 ] 


1.6 

1.5 

[g/m 2 ] 









*The material properties are given with reference to the ply coordinate axes where index 1 1 denotes the ply principal 
axis that coincides with the direction of maximum in-plane Young’s modulus (fiber direction). Index 22 denotes the 
direction transverse to the fiber in the plane of the lamina and index 33 the direction perpendicular to the plane of the 
lamina. 
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APPENDIX - B: FINITE ELEMENT MESH GENERATION 

Jointmesh - A Fortran based routine to generate Abaqus/Standard® input decks for 
configurations with overlaps and gaps 

A FORTRAN routine (jointmesh . f) was developed, which allows the automatic creation of an 
Abaqus/Standard® [24] input deck (.inp) and thus simplifies quick modifications in specimen 
dimensions and mesh size. The specimen dimensions and prescribed displacements as well as 
material data, ply orientation, and stacking sequence are provided by the user in a simple data file 
that is read by the routine at the beginning of the execution. The user can input up to 100 sections of 
arbitrary length in the y-direction, as shown in Figure Bl, where an example with 15 sections in the 
y-direction is shown. Please note that in the coordinate system used, x is the through-thickness 
direction and y is in length direction. For a three-dimensional problem, z is the width direction. The 
number of sections in the x and z-direction currently remain hard-coded and limited to three, 
however, the number of elements in all sections may be selected freely by the user. Thus, a 
multitude of zones along the length (y-direction) and through the thickness (x-direction) may be 
defined as shown in Figures Bl and B2. The user can assign arbitrary material properties (up to 
100) and orientations (up to 9) to each of the zones. The zones extend across the width (z-direction) 
as shown in Figure B2 for a three-dimensional model. The functions of the main program and the 
subroutines are briefly described in the following paragraphs together with the data input required 
from the user. An example data input file (mesh_input .data) and a shortened 
Abaqus/Standard® input deck (s_multi_gaps_linear . inp) are provided in the following 
sections. 


Main Program (JOINTMESH) 

The main program reads the input data from the file mesh_input . data, checks the validity of 
this data and makes the calls to the appropriate subroutines. As output, a file is created - 
qus/Standard® format which contains the nodal point coordinates, the element connectivity, material 
data as well as load and boundary condtions. The file can be used as input directly for a finite 
element analysis in Abaqus/Standard® [24], 

Data Input 

The data input file mesh_input . data is the only user required input. All variables which are 
necessary to generate the finite element model including the external loading, the boundary 
conditions and the material data are defined here. 

• nwrite : switch to obtain additional informational output to the screen during model 
generation 

• htext : problem name header, max. 60 characters 

• utext : unit header, max. 60 characters 

• abtext : name of output file to which model data will be written, max. 20 characters 

• el type : type of element used in the FE-model, max. eight characters 

• two dimensional four-node element in Abaqus/Standard® [24] 

• Abaqus/Standard® (plane-strain): CPE4, CPE4I, CPE4H, CPE4R, CPE4RH, 
CPEIH 

• Abaqus/Standard® (plane-stress): CPS4, CPS4I, CPS4R 

• two dimensional eight-node element 
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• Abaqus/Standard® (plane-strain): CPE8, CPE8R, CPE8H, CPE8RH 

• Abaqus/Standard® (plane-stress): CPS8, CPS8R 

• three dimensional eight-node brick element 

• Abaqus/Standard®: C3D8, C3D8H, C3D8I, C3D8IH, C3D8R, C3D8RH 

• three dimensional twenty-node element 

• Abaqus/Standard®: C3D20, C3D20H, C3D20R, C3D20RH 

• n : number of y- segments as shown in Figures B1 (15 here - max. 100) 

• YC (N) : array defining length (y) segments of specimen as shown in Figure B1 

• NY (N) : array defining number of elements in corresponding length segment 

• XB(3): array defining thickness (x) segments of specimen as shown in Figure B1 

• NX (3) : array defining number of elements in corresponding thickness segment 

• ZA(3) : array defining width (z) segments of specimen as shown in Figure B2 

• NZ (3) : array defining number of elements in corresponding width segment 

• nxyst(3) : array defining the start values for element zones 

• nxyend ( 3 ) : array defining the end values for element zones 

• z s t ( 3 ) : array defining start values for z 

• yst (3) : array defining start values for y 

• xs t ( 3 ) : array defining start values for x 

• with the input up to this point, the mesh is generated and the element sets in Abaqus/Standard® 
are defined as shown in Figure B2b. 

• matnum : number of materials in the model (a maximum of 100 materials can be 
defined using the following lines) 

• the following block of input has to be repeated matnum times (once for each material defined 
in the model) 

• mtext : material name. If the user chooses a value matnum>100 to act as a switch, a 
single line with the name of a file needs to be provided here which will be included in 
the Abaqus/Standard input file abtext. This file has to be a Abaqus/Standard input file 
. inp and contain all the necessary data to define a material, orientations and element 
sections. 

• naniso : switch to determine material type 

naniso=0 : isotropic material, two values E and v expected (ianiso=2) 

naniso=l : 2D laminate properties, four values expected (ianiso=4) 

naniso=2 : engineering constants, nine values expected (ianiso=9) 

naniso=3 : orthotropic material, nine values expected (ianiso=9) 

naniso=4 : anisotropic material, 21 values expected (ianiso=21) 

• matdat (ianiso) : material data, max. 21 values according to Abaqus/Standard® 

[24], 

• alpha (3): coefficients of thermal expansion for the material 

• be ta ( 3 ) : coefficients of moisture expansion for the material 

• no : number of orientations for the material (maximum 9) 

• the following block of input is repeated for each material 

• sys(6): orientation of the material's x' and y' direction 
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• nrot : axis about which the material layers are rotated 

• orient (no) : array containing angles of rotation defined for each material layer, 

max. 9 values 

• This ends the input for materials. If the user chooses a value matnum>100 to act as a switch, a 
single line with the name of a file needs to be provided here which will be included in the 
Abaqus/Standard® input file abtext. This file has to be a Abaqus/Standard® input file . inp 
and contain all the necessary data to define a material, orientations and element sections. 

• the following line of input has to be repeated nelx=nx (1) +nx (2) +nx (3) -times for each 
element layer in the x-direction 

• mzone (n) : array/line with n- values containing the material number and orientation 
number for each material zone in the model (e.g. 41 is material 4 with orientation 1; 10 
is material 1 with orientation 0) 

• thickn : thickness value for 2D type elements 

• nonlin : switch to perform non-linear analysis in Abaqus/Standard® for nonlin=l 

• i tnum : Abaqus/Standard®: maximum number of iterations 

• r(4) : Abaqus/Standard®: increment values 

• nhtl (3) : switch to activate analysis for mechanical, temperature and moisture loads 

• htl (5) : hygrothermal and mechanical load input 

■ initial temperature condition (Abaqus/Standard® input shown in Figure B3b) 

■ current temperature (Abaqus/Standard® input shown in Figure B3b) 

■ initial moisture condition (Abaqus/Standard® input shown in Figure B3c) 

■ current moisture condition (Abaqus/Standard® input shown in Figure B3c) 

■ applied displacement (Abaqus/Standard® input shown in Figure B3a) 

An example file is provided in the next section. 

Subroutine abelem 

This subroutine is used to convert the element names used in Abaqus/Standard® to internally used 
generic type numbers. 

call abelem (el type , nelpt , midnode , ntype , nwrite) 

• Input: 

• el type : element name or type in finite element code (maximum eight characters) 

• nwrite : switch to obtain additional informational output to the screen during model 
generation 

• Output: 

• ntype : internal element number 

• 204: two dimensional four-node element 

• 208: two dimensional eight-node element 

• 308: three dimensional eight-node element 

• 320: three dimensional twenty-node element 

• nelpt : number of nodes in this element as stored by the finite element code 

• midnode : for elements containing midside nodes midnode=2 
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Subroutine header 

This subroutine writes the header line to the Abaqus/Standard® input file (abtext). 

call header (ifecode ,htext,utext , iunit) 

• Input 

• ifecode : switch to determine which FE code will be used (this option is currently 
not used; if ecode=l hard-coded for Abaqus/Standard®) 

• htext : problem name header, max. 60 characters 

• utext : unit header, max. 60 characters 

• iunit : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 


Subroutine genpco 

This subroutine computes the nodal point coordinates for a three dimensional rectangular mesh in 
space. The *NODE keyword and the nodal point coordinates are written to the Abaqus/Standard® 
input file (abtext). 

call genpco (XB , YC , NX , NY , M , XST , YST , IUNIT , N , ZA , 

NZ , ZST , midnode , ifecode , nwrite ,ntype) 




Input 


XB (3) : array defining thickness (x) segments of specimen 

YC (5) : array defining length (y) segments of specimen 

NX (3) : array defining number of elements in corresponding thickness segment 

NY (5) : array defining number of elements in corresponding length segment 

M: size of array XB (M) , XST (M) ,YST(M) ,NX(M) 

XST (3) : start value 
YST (3) : start value 

IUNIT : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 

N : size of array yc(n), ny(n) 

ZA(3) : array defining width (z) segments of specimen 

NZ (3) : array defining number of elements in corresponding width segment 

ZST (3) : start value 

midnode : for elements containing midside nodes midnode=2 

ifecode : switch to determine which FE code will be used (this option is 

currently not used; if ecode=l hard-coded for Abaqus/Standard®) 
nwrite : switch to obtain additional output 
ntype : internal element number 
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Subroutine net2d 

This subroutine creates the finite element mesh by assigning nodal point numbers to a rectangular 
mesh consisting of four or eight-node 2D elements. The *ELEMENT keyword and the element 
connectivity are written to the Abaqus/Standard® input file (abtext). 


call ne t2d (NX , NY , M , N , IUNIT , ELNOD , NELPT , nxy s t , nxyend, midnode, 
ntype , eltype , if ecode , nwrite) ) 

• Input 

• NX (3) : array defining number of elements in corresponding thickness segment 

• NY (5) : array defining number of elements in corresponding length segment 

• M : size of array xb(m), xst(m),yst(m),nx(m) 

• N : size of array yc(n), ny(n) 

• IUNIT : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 

• ELNOD : array of element topology 

• NELPT : number of nodal points in element, size of array elnod 

• nxys t ( 3 ) : array defining the start values for element zones 

• nxyend ( 3 ) : array defining the end values for element zones 

• midnode : for elements containing midside nodes midnode=2 

• ntype : internal element number 

• eltype: type of element used in the FE-model, max. Eight characters 

• if ecode : if ecode : switch to determine which FE code will be used (this 

option is currently not used; if ecode=l hard-coded for Abaqus/Standard®) 

• nwrite : switch to obtain additional output 


Subroutine net3d 

This subroutine creates the finite element mesh by assigning nodal point numbers to a rectangular 
mesh consisting of eight or twenty-node 3D brick elements. The *ELEMENT keyword and the 
element connectivity are written to the Abaqus/Standard® input file (abtext). 

call ne t3d (NX , NY , NZ , M , N , IUNIT , ELNOD , NELPT , nxy s t , nxyend 
midnode , ntype , eltype , ifecode , nwrite) 

• Input 

• NX (3) : array defining number of elements in corresponding thickness segment 

• NY (5) : array defining number of elements in corresponding length segment 

• NZ (3) : array defining number of elements in corresponding width segment 

• M : size of array xb(m), xst(m),yst(m),nx(m) 

• N : size of array yc(n), ny(n) 

• IUNIT : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 

• ELNOD : array of element topology 
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• NELPT : number of nodal points in element, size of array elnod 

• nxys t ( 3 ) : array defining the start values for element zones 

• nxyend ( 3 ) : array defining the end values for element zones 

• midnode : for elements containing midside nodes midnode=2 

• ntype : internal element number 

• el type: type of element used in the FE-model, max. Eight characters 

• ifecode: ifecode: switch to determine which FE code will be used (this 

option is currently not used; if ecode=l hard-coded for Abaqus/Standard®) 

• nwrite : switch to obtain additional output 


Subroutine material 

This subroutine writes the material data *MATERIAL , *ELASTIC ... or 
*INCLUDE=matfile keywords to the Abaqus/Standard® input file (abtext). The 
♦ORIENTATION as well as *SOLID SECTION keywords are also written to the 
Abaqus/Standard® input file. 

call material (naniso , matdat , sys , nrot , orient , mm , nx , m , thickn , 
iuni t , ntype , matf ile , alpha , beta , mx , matnum , no , mzone , mtext , nn , N) 


Input 

• naniso : switch to determine material type 

• matdat : material data 

• sys (6): orientation of the material's x' and y' direction 

• nrot : axis about which the material layers are rotated 

• or ien t ( 2 0 ) : array containing the angles of rotation of each material layer, max. 

20 values 

• mm : number of layers or orientations 

• nx (3) : array defining number of elements in corresponding thickness segment 

• m : size of array xb(m), xst(m),yst(m),nx(m) 

• thickn : element thickness for 2D analysis 

• iunit : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 

• ntype : internal element number 

• matf ile : file name where Abaqus/Standard® material data are stored (for naniso=4), 
max. 20 characters 

• alpha (3) : coefficients of thermal expansion 

• beta (3) : coefficients of moisture expansion 

• mx : length of material array 

• matnum : number of materials used 

• no : number orientations defined for a material 

• mzone (n) : array/line with n- values containing the material number and 

orientation number for each material zone in the model 
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• mtext: material name 

• nn : maximum size of array yc(n), ny(n). Hard-coded n=100 

• N : actual size of array yc(n), ny(n) 

Subroutine abhist 

This subroutine sets up the Abaqus/Standard® “history” part of the input file (abtext). 

The* STEP, *END STEP as well as * STATIC keywords are written to the input file together 
with keywords that specify the hygro-thermal and mechanical loads. Keywords that specify the 
requested output are also written to the input file. 

call abhist (NX , NY , NZ , ZA , M , N , R , NFRD , IUNIT , nonlin , htl , 
nodef , i tnum , ntype , midnode , if ecode , nwri te , nhtl ) 

• Input 

• NX (3) : array defining number of elements in corresponding thickness segment 

• NY (5) : array defining number of elements in corresponding length segment 

• NZ (3) : array defining number of elements in corresponding width segment 

• ZA(3) : array defining width (z) segments of specimen 

• M : size of array xb(m), xst(m),yst(m),nx(m) 

• N : size of array yc(n), ny(n) 

• R(4) : Abaqus/Standard®: increment values 

• NFRD : number of degrees of freedom 

• IUNIT : fortran unit number to be supplied by the user (this option is currently not 

used; iunit=50 hard-coded). 

• nonlin : switch to non-linear analysis for nonlin=l 

• htl (5) : hygrothermal and mechanical load input 

• nodef : array 

• i tnum : Abaqus/Standard®: maximum number of iterations 

• ntype : internal element number 

• midnode : for elements containing midside nodes midnode=2 

• if ecode : switch to determine which FE code will be used 

• nwrite : switch to obtain additional output 

• nhtl (3) : switch to activate analysis for mechanical, temperature and moisture loads . 
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Example user input file: mesh_input . data 

This example file refers to the problem shown in Figure Bl. The input shown is for generating a 
model that contains 1 5 y-segments ( N= 1 5 ) and a total of 28 element layers (NX ( 3 ) = ( 2 , 24, 
2 ) ) across the thickness (x-direction). Four materials are defined for this example ( MATNUM=4 ) . 
Three material definitions are used to define the 0°, ±45° and 90° orientations of the IM-7/8552 tape 
and one material for the 8552 resin. For the definition of the IM-7/8552 tape plies, laminate 
properties are supplied (NANISO=l ) . The 8552 resin material is isotropic (NANISO=0 ) . Only 
for the IM-7/8552 tape plies, one orientation each has been defined (no=l ) . For each of the 28 
element layers across the thickness (x-direction), input is provided to assign a material and 
orientation to each of the y-segements along the length. The first digit represents the material and 
the second digit the orientation. Therefore, 31 denotes that material data from material 3 
(IM7/8552) and orientation 1 (±45°) are assigned to that particular section. Since material 4 (8552 
resin) is isotropic, no orientations have been assigned, and the sequence 4 0 is used. 


' NWRITE : WRITE SWITCH FOR OUTPUT COMMENT' 
’HEADER-TEXT: max 60 characters' 

'UNIT-TEXT: max 60 characters' 

'name of output file max 20 characters' 
'Element Type C3D20,CPE4, etc' 

'N: number of Y-SEGMENTS' 

'YC: LENGTH of Y-SEGMENTS' 

2.54 2.54 3.81 2.54 6.35 2.54 14.605 

'NY: NUMBER OF ELEMENTS IN Y-SEGMENTS' 

63 25 146 

'XB: length of x- segments' 

'NX: NUMBER OF ELEMENTS IN X- SEGMENTS' 

'ZA: LENGTH of z- segments' 

'NZ: NUMBER OF ELEMENTS IN Z - SEGMENTS ' 
'NXYST: first y-segment with elements' 
'NXYEND: last y-segment with elements' 

' ZST : start value ' 

' YST : start value' 

' XST : start value' 

'MATNUM: number of materials in the model' 

' material name ' 

' NANISO : =0 : iso , =1 : lam , =2eng , =3 : ort , =4 : ani ' 
'Material Input Ell, ... max. 21 values' 
'alpha: coefficients of thermal expansion' 
'beta: coefficients of moisture expansion' 

' number of orientations ' 

'SYS: material system, 6 values' 

'nrot: axis about which material is rotated' 
'orient: orientation of material in layers' 

' material name ' 

'NANISO: =0 : iso,=l : lam,=2eng,=3 : ort ,=4 : ani ' 
'Material Input Ell, ... max. 21 values' 
'alpha: coefficients of thermal expansion' 
'beta: coefficients of moisture expansion' 

' number of orientations ' 

'SYS: material system, 6 values' 

'nrot: axis about which material is rotated' 
'orient: orientation of material in layers' 

' material name ' 

'NANISO: =0 : iso,=l : lam,=2eng,=3 : ort ,=4 : ani ' 
'Material Input Ell, ... max. 21 values' 
'alpha: coefficients of thermal expansion' 
'beta: coefficients of moisture expansion' 

' number of orientations ' 

'SYS: material system, 6 values' 

'nrot: axis about which material is rotated' 


3 

'section - delta=6.35 mm, c=2.54 mm, multi gaps' 
'units: mm, N, MPa' 

' s_multi_gaps_l inear . inp ' 

' CPS8 ' 

15 

14.605 2.54 6.35 2.54 3.81 2.54 2.54 6.35 

146 25 63 25 144 96 96 240 96 96 144 25 

0.3 0.6 0.3 

2 24 2 

5.0 15.4 5.0 

10 15 10 

111 
15 15 15 

0 . 0 . 0 . 

76.2 0. 0. 

0 . 0 . 0 . 

4 

' IM7/8552-0 ' 

1 

160523. 14386. 0.467 5171. 

-9.E-7 2.88E-5 2.88E-5 


0. 

3310. E-6 

3310. E-6 

1 

0.0 

o 

o 

o 

1—1 

1 

1.0 0.0 


3 

0. 

' IM7/8552-90 ' 

1 

9809. 12422. 0.461 3923. 

2.88E-5 2.88E-5 -9.E-7 

3310. E-6 3310. E-6 0. 

1 

0.0 - 1.0 0.0 1.0 0.0 0.0 
3 
0 . 

' IM7/8552-45 ' 

1 

32381. 13810. 0.464 4547. 

1.39E-5 2.88E-5 1.39E-5 

1655. E-6 3310. E-6 1655. E-6 

1 

0.0 - 1.0 0.0 1.0 0.0 0.0 
3 
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' orient : orientation of material in layers ' 
'material name' 

' NANISO : =0 : iso , =1 : lam , =2eng , =3 : ort , =4 : ani ' 
'Material Input Ell, ... max. 21 values' 
'alpha: coefficients of thermal expansion' 
'beta: coefficients of moisture expansion' 

' number of orientations ' 

'material zones x(l) ' 

'material zones x(2) ' 

'material zones x(3) ' 

'material zones x(4) ' 

'material zones x(5) ' 

'material zones x(6) ' 

'material zones x(7)' 

'material zones x(8) ' 

'material zones x(9) ' 

'material zones x(10)' 

'material zones x(ll)' 

'material zones x(12)' 

'material zones x(13) ' 

'material zones x(14)' 

'material zones x(15)' 

'material zones x(16)' 

'material zones x(17)' 

'material zones x(18)' 

'material zones x(19) ' 

'material zones x(20)' 

'material zones x(21)' 

'material zones x(22)' 

'material zones x(23) ' 

'material zones x(24)' 

'material zones x(25)' 

'material zones x(26)' 

'material zones x(27)' 

'material zones x(28)' 

' thickness for 2d elements , specimen width ' 
'NONLIN: =1 for non-linear analysis' 

'itnum: maximum numer of iterations' 

'R: 4 increment values' 

'nhtl: mechanical, temp, hygro - on (1) ,off (0) ' 
'htl: inti cond, temp, inti cond, hygro, disp' 


0. 

' R8552 ' 

0 

11380. 0.3 
2.88E-5 
3310. E-6 

0 

31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 21 21 21 21 21 21 21 21 21 21 40 21 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

31 31 31 31 31 31 31 31 31 31 31 40 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

11 11 11 11 11 11 11 11 40 31 31 31 31 31 31 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

31 31 31 31 31 31 40 11 11 11 11 11 11 11 11 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

21 21 21 40 31 31 31 31 31 31 31 31 31 31 31 

31 40 21 21 21 21 21 21 21 21 21 21 21 21 21 

31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 

25.4 
0 

500 

0.1 1 . 0.001 1 . 

10 0 

450. 293. 0. 1.0 0.762 


Example Abaqus/Standard® input file: s_multi_gaps_linear . inp 

The FORTRAN routine Jointmesh discussed above generates Abaqus/Standard® input file (.inp) 
which can be used directly to perform a finite element analysis in Abaqus/Standard®. It is assumed 
that the reader is familiar with Abaqus/Standard® and the syntax used in the input file. The 
example input file is given to provide an overview of an entire analysis and assist the readers in 
creating their own analyses. 

♦HEADING 

section - delta=6.35 mm, c=2.54 mm, multi gaps 
units: mm, N, MPa 
*NODE , NSET=NALL 

1 , 0.0000000E+00 , 76.20000 , 0 . 0000000E+00 

2 , 7 . 5000003E-02 , 76.20000 , 0 . 0000000E+00 

3 , 0.1500000 , 76.20000 , 0 . 0000000E+00 

*NSET , NSET=NFORCE , GENERATE 

1 , 57 , 1 

*NSET , NSET=NBOUNDL , GENERATE 

1 , 57 , 1 

*NSET , NSET=NBOUNDR , GENERATE 

163021 , 163077 , 1 

*ELEMENT , TYPE=CPS8 , ELSET=EALL 

1 , 1 , 115 , 

58 , 116 , 60 , 
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*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 
*ELSET , 

*ELSET , 


2 

172 


115 , 229 , 231 , 

230 , 174 , 116 


ELSET=LY01, GENERATE 

1 , 1430 , 1 

ELSET=LY01-01, GENERATE 

1 , 146 , 1 

ELSET=LY01-02 , GENERATE 
147 , 171 , 1 

ELSET=LY01-03 , GENERATE 
172 , 234 , 1 

ELSET=LY01-04 , GENERATE 
235 , 259 , 1 

ELSET=LY01-05 , GENERATE 
260 , 403 , 1 

ELSET=LY01-06, GENERATE 
404 , 499 , 1 

ELSET=LY01-07 , GENERATE 
500 , 595 , 1 

ELSET=LY01-08 , GENERATE 
596 , 835 , 1 

ELSET=LY01-09, GENERATE 
836 , 931 , 1 

ELSET=LY01-10 , GENERATE 
932 , 1027 , 1 

ELSET=LY01-11 , GENERATE 
1028 , 1171 , 1 

ELSET=LY01-12 , GENERATE 
1172 , 1196 , 1 

ELSET=LY01-13 , GENERATE 
1197 , 1259 , 1 

ELSET=LY01-14 , GENERATE 
1260 , 1284 , 1 

ELSET=LY01-15, GENERATE 
1285 , 1430 , 1 

ELSET=LY02 , GENERATE 
1431 , 2860 , 1 

ELSET=LY02-01, GENERATE 
1431 , 1576 , 1 

ELSET=LY02-02 , GENERATE 
1577 , 1601 , 1 

ELSET=LY02-03 , GENERATE 
1602 , 1664 , 1 

ELSET=LY02-04 , GENERATE 
1665 , 1689 , 1 

ELSET=LY02-05 , GENERATE 
1690 , 1833 , 1 

ELSET=LY02-06 , GENERATE 
1834 , 1929 , 1 

ELSET=LY02-07 , GENERATE 
1930 , 2025 , 1 

ELSET=LY02-08 , GENERATE 
2026 , 2265 , 1 

ELSET=LY02-09, GENERATE 
2266 , 2361 , 1 

ELSET=LY02-10 , GENERATE 
2362 , 2457 , 1 

ELSET=LY02-11, GENERATE 
2458 , 2601 , 1 

ELSET=LY02-12 , GENERATE 
2602 , 2626 , 1 

ELSET=LY02-13 , GENERATE 
2627 , 2689 , 1 

ELSET=LY02-14 , GENERATE 
2690 , 2714 , 1 

ELSET=LY02-15 , GENERATE 
2715 , 2860 , 1 


ELSET=LY28-15, GENERATE 


117 , 
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40040 


1 


39895 , 

*BOUNDARY 
NBOUNDR, 1, 2 
NBOUNDL , 1, 1 
*Initial Conditions, type=TEMPERATURE 
NALL, 450.00000 

*Initial conditions, type=field, variable=l 
NALL, 0.00000 

*MATERIAL , NAME=IM7/ 8552-0 
*ELASTIC , TYPE=LAMINA 

0 . 16052E+06 , 14386. , 0.46700 , 5171.0 

*Expansion, type=ortho 

-8 . 9999997E-07 , 2 . 8799999E-05 , 2 . 8799999E-05 

*Expansion, type=ortho, field=l 

0.0000000E+00 , 3 . 3100001E-03 , 3 . 3100001E-03 

* ORIENTATION , NAME=LY0011, SYSTEM=RECTANGULAR , DEFINITION=COORDINATES 

0.0000 ,- 1.00000 , 0.0000 , 1.00000 , 0.0000 , 0.0000 

3, 0.000 

*MATERIAL , NAME=IM7/8552-90 
*ELASTIC , TYPE=LAMINA 

9809.0 , 12422. , 0.46100 , 3923.0 


*Expansion, type=ortho 

2 . 8799999E-05 , 2 . 8799999E-05 , -8 . 9999997E-07 

*Expansion, type=ortho, field=l 

3 . 3100001E-03 , 3 . 3100001E-03 , 0 . 0000000E+00 

* ORIENTATION , NAME=LY0021, SYSTEM=RECTANGULAR , DEFINITION=COORDINATES 

0.0000 ,- 1.00000 , 0.0000 , 1.00000 , 0.0000 , 0.0000 

3, 0.000 

*MATERIAL , NAME=IM7/8552-45 
*ELASTIC , TYPE=LAMINA 

32381. , 13810. , 0.46400 , 4547.0 


*Expansion, type=ortho 

1 . 3900000E-05 , 2 . 8799999E-05 , 1 . 3900000E-05 

*Expansion, type=ortho, field=l 

1.6550000E-03 , 3 . 3100001E-03 , 1 . 6550000E-03 

* ORIENTATION , NAME=LY0031, SYSTEM=RECTANGULAR , DEFINITION=COORDINATES 
0.0000 ,- 1.00000 , 0.0000 , 1.00000 , 0.0000 , 0.0000 
3, 0.000 

*MATERIAL , NAME=R8552 
*ELASTIC , TYPE=ISOTROPIC 

11380.00 , 0.3000000 

*Expansion, type=iso 
2 . 8799999E-05 


*Expansion, type=iso, field=l 
3 . 3100001E-03 



*SOLID SECTION, 
25.40000 

ELSET=LY01-01, 

ORIENTATION=LY0031 , 

MATERIAL=IM7/8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-02, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-03, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-04, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-05, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY0 1 -06, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-07, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY0 1 -08, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY0 1 -09, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-10, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY0 1 -11, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-12, 

ORI ENTAT I ON=L Y 0031 , 

MATERI AL= IM7 /8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-13, 

ORIENTATION=LY0031 , 

MATERI AL= IM7 /8552-45 
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*SOLID SECTION, 
25.40000 

ELSET=LY01-14, 

ORI ENTAT I ON=L Y 0031 , 

MATERIAL=IM7/8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY01-15, 

ORI ENTAT I ON=L Y 0031 , 

MATERIAL=IM7/8552-45 

*SOLID SECTION, 
25.40000 

ELSET=LY02-01, 

ORI ENTAT I ON=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-02, 

ORI ENTAT I ON=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-03, 

ORI ENTAT I ON=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-04, 

ORI ENTAT I ON=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-05, 

ORIENTAT ION=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-06, 

ORIENTAT ION=L Y 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-07, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-08, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-09, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-10, 

ORIENTAT ION=LY 0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-11, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-12, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-13, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY02-14, 

MATERI AL=R8 552 


*SOLID SECTION, 
25.40000 

ELSET=LY02-15, 

ORIENTAT ION=LY0021 , 

MATERIAL=IM7/8552-90 

*SOLID SECTION, 
25.40000 

ELSET=LY2 8-15, 

ORIENTAT ION=LY00 31 , 

MATERIAL=IM7/8552-45 


★STEP 

*STATIC 

1 . 0000 , 1 . 0000 , 

*BOUNDARY 

NFORCE , 2, 2, 0.76200 

*OUTPUT , FIELD, VARIABLE=PRE SELECT , FREQ=1 
*ELEMENT OUTPUT, ELSET=EALL 
TEMP, FV1 

* Output , history, VARIABLE=PRESELECT , f req=l 

*NODE output , NSET=NBOUNDR 

RF2 

*NODE output , NSET=NFORCE 
RF2 

*NODE output , NSET=NBOUNDR 
U2 

*NODE output , NSET=NFORCE 
U2 

*END STEP 
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Figure B1 . Example of input to user datafile mesh_input.data. 



a. Finite element mesh. 



b. Detail of local zones. 


Figure B2. Example of a three-dimensional model with element sets to define 

material properties. 
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x,u,1 



* BOUNDARY 
NBOUNDR , 1, 2 
NBOUNDL ,1,1 


* BOUNDARY 

NFORCE , 2, 2, 0.10000 


U-COnst *NSET , NSET-NFORCE *NSET NSET=NBOUNDR U_ ® 

V=0 *NSET , NSET=NBOUNDL V=0 

a. Boundary conditions and mechanical load at node sets. 



A 

* 

A 

$ 

A 

$ 

A 

$ 

A 

¥ 


TEMP 
(Avg: 75%) 

B +5.000e+02 
+4.811e+02 
+4.622e+02 
+4.432e+02 
— I- +4.243e+02 

+4.054e+02 

+3.865e+02 

+3.676e+02 

+3.487e+02 

+3.297e+02 

H +3.1 086+02 
+2.919e+02 
+2.730e+02 


♦Initial Conditions, type=TEMPERATURE 
NALL, 450.00000 

♦Temperature 
NALL, 293.00000 


x,u,1 


AT=157K=const. 

difference between cure and room temperature 
b. Thermal load defined as field. 


y,v,2 


FV1 

(Avg: 75%) 

B +1 ,200e-02 
+1 . 1 00e-02 
+1 .000e-02 
+9.000e-03 
— I- +8.000e-03 

+7.000e-03 

+6.000e-03 

+5.000e-03 

+4.000e-03 

+3.000e— 03 

H +2.0006—03 
+1 .000e-03 
+0.000e+00 


*Initial conditions, type=field, variable=l 
NALL, 0.00000 

*field, variable=l 
NALL, 0.01000 


X,U,1 


c=0.01 (1% mass/weight)=saturation y,v,2 

moisture diffusion through the thickness of the 
specimen not computed at this point 

c. Moisture load defined as field. 

Figure B3. Abaqus input generated for load and boundary conditions. 
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